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This  note  deals  with  the  integral 
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.2  J (rx)Y  (ax)-J  (ax)Y  (rx)  dx 
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which  arises  in  the  solution  of  the  equation  of  heat  flow 
in  a region  internally  bounded  by  a circular  cylinder  with 


constant  surface  temperature  [1,  pp.  280-282],  Its  eval- 
uation for  large  values  of  t presents  certain  difficulties 
because  the  standard  Abelian  theorems  for  Laplace  integrals 
are  not  applicable.  Although  asymptotic  expansions  of 
similar  integrals  can  be  found  in  the  literature  [1,3?*+]? 
they  generally  are  stated  without  explanation  of  their  de- 
rivation, It  seems  desirable,  therefore,  to  bridge  this 
gap  by  giving  a rather  complete  discussion,  which  may 
prove  useful  in  the  evaluation  of  related  integrals. 

The  problem  was  brought  to  our  attention  by  B.  Peavy  of 
the  National  Bureau  of  Standards;  we  appreciate  his  interest 
in  the  publication  of  our  result  in  this  form. 


1 , Let 


*This  paper  was  prepared  under  a National  Bureau  of  Stand- 
ards contract  with  The  American  University. 
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(1  ) F(x)  = 


and  put 
(2)  f(t) 


J (rx)Y  (ax)-J  (ax)Y  (rx) 
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[ Jo2(ax)+Yo2(ax) ] ,x 
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= JT  r e kx  tF(x)dx 


where  it  is  assumed  that  a,  k,  r are  positive  constants 

From  the  definitions  of  the  Bessel  functions  J (z), 

o 

Y (z)  we  obtain  the  expressions 


(3) 


(4) 


J (rx)Y  (ax)  - J (ax)Y  (rx)  = “ ~ [ 1 + 0(x2)], 

O O O O **  L 

J 2(ax)  + Y 2 (ax)  = [1  + (~  £ n ^9ax)2][  1 + 0(x2)]* 


where,  as  customary,  0(x  ) denotes  a function  such 

tha  t|pb(x) /x2  < A for  some  A > 0 as  x — ^0.  Consequently, 

we  find 


(5) 


F (x)  = % in  f [1  + (2/„  -to)2]-1  ) 1 = 

r 1 + O(x^)  XL 

£ n “ ([1t(/n  “^-7p)2]x)  1 [ 1 + G(x)  ] 
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where  the  function  G(x)  is  such  that  G(x)  = 0(x  ) as  x — ?0 
and  G(x)  - 0(ex)  as  x — ?oo#  Thus,  F(x)  has  an  infinite 
discontinuity  at  x = 0,  and  the  integral  in  (2)  is  for 
any(finite)  value  of  t an  improper  integral.  However,  if 


*Throughout , C~=~ oT577»  « . Euler’s  constant 
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we  denote  by  g(x)  the  numerator  in  the  second  expression  of 
(5)  and  note  that  |g(x)|  <_ M on  any  interval  0 <_ x <_X,  then 
for  any  sufficiently  small  > 0 


so  that  the  integral  in  (2)  certainly  converges  at  its  lower 
limit  for  any  value  of  t _>  0.  The  convergence  at  its  upper 
limit  is  evident. 

It  is  to  be  noted  that  the  fact  that  the  integral  in  (2) 
is  an  improper  integral  a fortiori  precludes  the  possibility 
of  using  directly  numerical  methods  for  its  evaluation. 

And  this,  of  course,  holds  true  regardless  of  the  value  of 
t for  which  (2)  is  to  be  calculated. 

2.  In  analogy  to  the  asymptotic  theory  of  functions 
defined  by  Laplace  integrals  it  can  be  shown  that  the  im- 
portant contribution  to  f(t)  for  large  values  of  t will 
come  from  the  values  of  F(x)  near  x = 0.  Hence  we  need  to 
consider  the  two  integrals,  obtained  by  substituting  the 
expression  (5)  into  (2),  namely. 
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= f 1 (t)  + f2(t) . 


We  shall  first  deal  with  the  second  integral  f2(t) 


2 . 


Putting  y = kx  in  f2(t)  we  find 


00 

(8)  f2(t)  = ^ ^ J e ytG(v'y7k)[l  +(~£n  ^2°  7|)2]  1 ^ . 

The  integral  on  the  right  is  clearly  a Laplace  integral  to 
which  we  can  apply  the  standard  Abelian  theorems  of  the 
asymptotic  theory  of  Laplace  integrals  [2,  p,  202],  We  con- 
clude that  for  t — >oo 

(9)  f2( f B[1  + (f/nVf)2]"1  (1/t) 

p 

where  c = [7^a/(2>/k)]  and  D is  the  constant  implied  by 
G(vyA)  = 0(y)  as  y — 7>0. 

3.  We  now  come  to  the  much  longer  discussion  of  the  in- 
tegral I(t)  in  f^(t),  i.e, 

(10)  I(t)  = J e“kx2t[l  + (fin  92]-1  f . 
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Putting  kx  = (1/c  )y  where  c = [?fa/(2>/k)]  we  find 

oo.  t 


(11 ) I(t)  = i I e 


e [1  + Ciiny)2]-1  f = 
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e c d (arc  tan  9 


Performing  an  integration  by  parts  on  the  last  integral,  we 

obtain 
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We  now  let  ty  = czj  noting  that  J'  e "'dz  = 1 , we  then  can 

o 

oo 


write 


7T*“  7T 


Kt)  = ip  + o 


e ''arctan  [ v^n(pz/t)  ]^2  _ 
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e ^dz  + P e zarctan[ ^ ^z 


oo 
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f f e_?[f  + arctan[^^])dz. 


Recalling  that  arctan  (oo)  = tt/2  and  making  use  of  an  elem- 


entary trigonometric  identity,  we  finally  obtain  from  (13) 
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(14) 


I ( t ) = — C e 7arctao  L "iT/ ”^)J  dz.* 

J o 

In  this  integral  we  pot  s = Lw  — which  yields 

OO 

(l£)  I(t)  = Q(s)  ~ ~~  2 ^ e 7!arctan  i^s/Cs/n?  - 1)j  dz 


TT 

2 


00 

e zarctan  [ ?rs/(  1 - s/n?)]  d: 


It  is  this  last  integral  which  we  shall  evaluate  for  small  s 
4.  To  obtain  an  asymptotic  representation  of  (15)  for 
s — ^0  we  shall  prove  that 


(16)  lim 
s— >0 


1 °o 
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e a re  tan  [ tr  s / ( 1 - s 


= irt 


Vi/e  must  show  that  given  any  > 0,  however  small,  there 
exists  a ^ (£)  > 0 such  that 


07) 


00 

f e 7arctan  [irs/(1  - $J!n z)]dz  - 7f  < £ 


for  all  s < d (£)  . Evidently,  this  inequality  will  be  sat- 
isfied if  for  all  s < ^ (£ ) 


*The  use  of  the  trigonometric  identity  involving  arctangents 
was  pointed  out  to  us  by  P.  ilenrici.  Our  original  argument 
wa s so raewh at  longer. 
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arctan  [^s/(1  - s^rz)]  - 'r'|  dz| 


(18)  0 
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| ^ e 7 “ arctan  [?rs/(1  - s^nzOdzj  + ir  | ^ e ?dzj  < £ 


where  7 > 0 is  some  constant  at  our  disposal.  Since  the 

function 


(19) 


h(t,z)  = arctan  Ort/O  - t^/nz)] 


is  uniformly  bounded  for  all  7 and  all  t > 0 and  h(t?z} — 
as  t — ^0  for  z finite,  we  have 

00  00 


(20)  j e 7 ^ arctan  [tts/(1  - s /n?)  ]dz 
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< K J e 7d7.  = Ke 


Choose  Z > 0 so  large  that  max  (Le  v,  we  <J)  <_  — 9 Then  all 

we  have  to  prove  is  that  there  exists  a ^ (£)  > 0 such  that 


(21) 


ll 


e z )“  arctan  l^s/O  - s 


/nz)  j - 7r  | dzj  <_ 


for  all  s (£) . This,  however,  is  now  trivial  since  the 

boundedness  of  h(t,z)  and  h ( t , z)  - — ?0  as  t — >0  together 
imply  the  bounded  convergence  of  (21 ) to  zero  as  s — * 0o 
By  virtue  of  (16)  we  can  conclude  from  0 5)  that  for 
s — * 0 
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(22) 


9(s) 


7 r 


whence  for  t — * oo 


(23) 


I(l)~  2~£MT/7) 


Hecalling  that  we  have  from  (7)  and  (10) 


(24) 


f1  (t)  = (~-2-  Jin  ~)  I(t) 


7T 


we  finally  obtain 
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>^n4t  k - 2^nTf  a 


If  we  compare  this  relation  with  that  derived  in  (9)  for 
f^t),  we  see  that  for  sufficiently  large  values  of  t the 
contribution  of  to  can  neglected. 

Thus,  the  behavior  of  f(t)  for  large  values  of  t is 
described  by  the  asymptotic  relation 


5.  It  is  possible  to  derive  similarly  higher  order 
terras  in  the  asymptotic  relation  (25)  for  f^(t).  In  fact 
we  can  obtain  in  this  manner  the  asymptotic  expansion 


for  t — *oo  . In  view  of  the  asymptotic  relation  (9)  for 
f2(f)  this  expansion  can  be  considered  as  an  asymptotic 


(26) 


f ^^n4tk  - 2^n~tfa 
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expansion,  in  the  sense  of  Poincare'",  of  the  function  f(t) 
for  t — >oo  , 

In  conclusion  we  remark  that  the  expansion 
oo 


(28) 
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J^(ax)+Y^(ax)  x 2 n*+tk  ~ 2/nU 
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(^n^tk  - 2^n  ^a) 


+ ♦ • » ] 


can  be  found  in  the  literature  [1,3AJ?  where  it  is  stated 
without  derivation. 


REFERENCES 


1.  Carslaw,  H.5.  and  Jaeger,  J.C.  Conduction  of  Heat  in 

Solid  s . Oxford  University  Pre¥s^  19^+7. 

2.  Doetsch,  G.  Theorie  und  Anwendunq  der  Laolace-Transf ornia- 

tion . (American  Ed.).  Dover  Publications,  1 9^+X. 

3.  Goldstein  S.  Some  two-dimensional  Diffusion  Problems  with 

Circular  Symmetry.  Proc.  London  Math.  Soc.,  vol. 

3^0932),  pp.  51-88. 

h.  Jaeger,  J.  C.  Heat  Flow  in  the  Pieoion  bounded  Internally 
by  a Circular  Cylinder.  Proc.  Roy.  Soc.  Edinburgh 
(A),  vol.  61 (1942),  pp.  223-228. 


The  American  University 


■ 


THE  NATIONAL  BUREAU  OF  STANDARDS 


Functions  and  Activities 

The  functions  of  the  National  Bureau  of  Standards  are  set  forth  in  the  Act  of  Congress,  March 
3,  1901,  as  amended  by  Congress  in  Public  Law  619,  1950.  These  include  the  development  and 
maintenance  of  the  national  standards  of  measurement  and  the  provision  of  means  and  methods 
for  making  measurements  consistent  with  these  standards;  the  determination  of  physical  constants 
and  properties  of  materials;  the  development  of  methods  and  instruments  for  testing  materials, 
devices,  and  structures;  advisory  services  to  Government  Agencies  on  scientific  and  technical 
problems;  invention  and  development  of  devices  to  serve  special  needs  of  the  Government;  and  the 
development  of  standard  practices,  codes,  and  specifications.  The  work  includes  basic  and  applied 
research,  development,  engineering,  instrumentation,  testing,  evaluation,  calibration  services,  and 
various  consultation  and  information  services.  A major  portion  of  the  Bureau’s  work  is  performed 
for  other  Government  Agencies,  particularly  the  Department  of  Defense  and  the  Atomic  Energy 
Commission.  The  scope  of  activities  is  suggested  by  the  listing  of  divisions  and  sections  on  the 
inside  of  the  front  cover. 

Reports  and  Publications 

The  results  of  the  Bureau’s  work  take  the  form  of  either  actual  equipment  and  devices  or 
published  papers  and  reports.  Reports  are  issued  to  the  sponsoring  agency  of  a particular  project 
or  program.  Published  papers  appear  either  in  the  Bureau’s  own  series  of  publications  or  in  the 
journals  of  professional  and  scientific  societies.  The  Bureau  itself  publishes  three  monthly  peri- 
odicals, available  from  the  Government  Printing  Office:  The  Journal  of  Research,  which  presents 
complete  papers  reporting  technical  investigations;  the  Technical  News  Bulletin,  which  presents 
summary  and  preliminary  reports  on  work  in  progress;  and  Basic  Radio  Propagation  Predictions, 
which  provides  data  for  determining  the  best  frequencies  to  use  for  radio  communications  throughout 
the  world.  There  are  also  five  series  of  nonperiodical  publications:  The  Applied  Mathematics 
Series,  Circulars,  Handbooks,  Building  Materials  and  Structures  Reports,  and  Miscellaneous 
Publications. 

Information  on  the  Bureau’s  publications  can  be  found  in  NBS  Circular  460,  Publications  of 
the  National  Bureau  of  Standards  ($1.00).  Information  on  calibration  services  and  fees  can  be 
found  in  NBS  Circular  483,  Testing  by  the  National  Bureau  of  Standards  (25  cents).  Both  are 
available  from  the  Government  Printing  Office.  Inquiries  regarding  the  Bureau’s  reports  and 
publications  should  be  addressed  to  the  Office  of  Scientific  Publications,  National  Bureau  of  Stand- 
ards, Washington  25,  D.  C. 


r 


